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Abstract 

The van Dam-Veltman-Zakharov (vDVZ) discontinuity requires that the mass m of 
the graviton is exactly zero, otherwise measurements of the deflection of starlight 
by the Sun and the precession of Mercury's perihelion would conflict with their 
theoretical values. This theoretical discontinuity is open to question for numerous 
reasons. In this paper we show from a phenomenological viewpoint that the m > 
hypothesis is in accord with Supernova la and CMB observations, and that the large 
scale structure of the universe suggests that m ~ 10~ 30 eV/c 2 . 
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1 The vDVZ discontinuity 



Recent papers have revisited the question of whether gravity could be medi- 



ated by an ultralight but massive graviton ( 


Will. 


2001 


Creminelli et al. 


2005; 


Deffavet & Rombouts. 


2005; 


Goldhaber & Nieto. 


2009; 


Arun & Will. 


2009). 


The classic works of 


van Dam & Veltman 


(1970) andi 


?akharov 


(1970) demon- 



strated that, in a linearized theory of massive gravity, the tensor gravitons of 
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General Relativity (GR), which couple to Tjk, the stress-energy tensor, are 
supplemented by a scalar graviton, which couples to Tjj, the trace of the stress- 
energy tensor. In the open limit as m — > 0+, the gravitational forces between 
non-relativistic masses can be accom odated simply by inc r easin g the gravita- 
tional constant G by the factor 4/3 (jMaggiord . 120081 ; IZed . 120031 ); whereas for 
m — 0, there is no scalar graviton hence G is not increased to compensate for 
the discontinuity An increased G for m > would lead to a prediction for the 
angle starlight is deflected by the Sun that is 4/3 that of GR. The agre ement 
betw een GR and this key observational test is now at the ~ 10 -5 level (jWilll . 
20011 ). which implies that the graviton must be strictly massless. 



Various elaborat e explanations for re solving the vDVZ disc o ntinuity have bee n 
suggested (e.g.. IVainshteinl fll972h : IPeffavet et ail <j2002h : iGruzinov! fj2005h : 
Gabadadze fc Gruzinovl (120051) ). In the present p a per w e follow the approach 
of iGoldhaber fc Nietol (l2009h and lArun fc Will d2Q09h by sidestepping the 
theoretical question of whether and how a theory of massive gravitons can be 
formulated, and focusing instead on observational data. Specifically we show 
that Supernova la and CMB data support the m > hypothesis, and that the 
large scale structure of the universe suggests that m ~ 10~ 30 eV/c 2 . 



2 A cosmological expansion discontinuity 



Consider a spatially uniform density (p = p(t)) expanding universe for which 
the spherical coordinate line element is of the form 



ds 2 = c 2 dt 2 - a 2 (t) [dr 2 + f 2 (r) dQ], 



From the null geodesic (ds = 0) for light traveling along a radial path (dQ = 0), 
we then have 



cdt = a dr. 



(2) 



The Einstein-de Sitter flPeeblesl . Il993h universe (GR with cosmological con- 
stant A = 0, f(r) — r, m — 0) requires that a oc t 2//3 , so if a photon is emitted 
at t e when a(t e ) = a e and observed at t when a(t„) = a„, the distance between 
the emitter (e) and observer (o) at time t is (jPeeblesl . 119931 ) 



a r = (2c/H )[l-(l + z) 



-1/21 



(3) 



where H Q = H(t a ) is Hubble's constant, and z = a /a e — 1 is the ob served 
redshift. If, however, m > 0, then Birkhoff's theorem (jPeeblesl . 119931 ) does 



2 



not apply, so we cannot model the gravitational field on the surface of a 
sphere without considering its surroundings, even if (as in this model) the 
surroundin gs have a constan t density. Instead, we adopt the phenomenological 
approach of lMaggior 2 (l2008h who hypothesized that the gravitational potential 



energy of two point masses, mi and m 2 , separated by the distance r is 
Gm 1 m 2 

V Y {r) = exp(-/xr), (4) 



where p = mc/h. Maggiore, following the IVainshteinl (119721 ) scenario, used 
G = (4/3)£tjv, where Gn is the Newtonian gravitational constant, but we 
shall ignore any distinction (if indeed one exists) between G and Gn- Then 
for a uniform density p universe, the gravitational potential at r = is formally 



Y 



-AnGp 



exp(— /ir) 



r 2 dr = —A-nGpX" 



-3GM A /A, 



(5) 



where X = 1/p is the reduced Compton wavelength of the graviton, and M\ 
is the mass of a sphere with radius X. (This outwardly naive solution does not 
apply if m = 0. See Appendix |A] for its rigorous derivation and a discussion 
of the m = discontinuity.) The integral of Eq. [5] is independent of where we 
select the origin, so the potential has no gradient and there is no cosmological 
acceleration; that is, a = constant and Eq. [3] is not valid. On a cosmological 
scale this is a Milne universe, for which the Einstein tensor G^ must be zero. 
To effect this, we set 

a{t) = a(0) + at, (6) 



and, in Eq. [fl 



f(r) = — sinh — . 
a c 



(7) 



One can easily verify that Eqs. ED [6] and 7\ defin e a metric for which G^ = 
by using Parker's Mathematica (Wolfram . 19991 ) notebook, Curvature and the 
Einstein Equation ( jHartld . 120021 ) . Appendix IB1 details how Eq. [7|was derived 
and verified. 



For m > 0, the integral of Eq. [2] is 



t„ 



cdt c f da c c . 

— - j — = — ln(a /a e ) = - m(l + z). 

a a J a a a 

a. 
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On substituting a = a H , we find that the distance between the emitter and 



observer at t is 



a a r = (c/H ) ln(l + z). 



(9) 



The area of the wavefront at t Q is 4n(a f) 2 , where 



c c 
a f = w sinh [ln(l + z)] = — - 

tl n & tl n 



1+2- 



1 + z 



(10) 



which, except for the sign, is in agreement with ISethi. Dev fc Jainl (120051 ) for 
a linear coasting cosmology (their Eq. 8, with Z\ — > z, z 2 — ► 0). 



3 The agreement with the type la supernovae (SNe la) observa- 
tions 



Two large research teams ( Riess et al. . 1998 ; Perlmutter et al. . 19991 ) recently 
found discrepancies between the distances to high redshift [z ~ 1) Type la 
supernovae (considered to be well understood "standard candles") when the 
SNe la distances are determined by their apparent magnitudes versus when 
they are determined by their redshifts using Eq. [3j The SNe la brightnesses 
appear to be about 25% weaker than expected, and so their distances are cor- 
respondingly greater than their redshifts would indicate. To phrase it another 
way, the SNe la redshifts are smaller than their magnitudes would indicate, 
and so a in the past appears to be smaller than it is today; that is, a > 
and the expansion rate of the universe is accelerating. Actually, however, the 
acceleration satisfies a > a mo dei oc — 1 -4//3 if the model is that of an Einstein- 
de Sitter universe. Although a is not necessarily greater than zero at present, 
there still is an unmodeled effect, an ostensibly repulsive force that has been 
ascribed to "dark energy" , that enters the GR field equations as a non-zero A. 
We offer an alternative to dark energy, namely that gravitons are not massless. 



In Fig . [J we compare the most recent set of high-confidence ("gold") (IRiess et al. 



20071 ) data for 182 SNe la sources distributed over the interval 0.0233 < 
z < 1.755 with two models in which A = 0: GR for a flat universe, with 
m = (Eq. [3]), and GR with m > (Eq. [9]). The Hubble constant was de- 
termined by a weighted least-squares adjustment of the model to the data. 
For m > 0, H Q = 60.1 km s _1 Mpc^x 2 = 179), and for m = 0, H = 
55.4 km s _1 Mpc _1 (% 2 = 288). Because the number of degrees of freedom 
v = 180 is large, x = (y 2 — u)/\/2u is approximately normall y distributed 



with unit variance (jCramerl . Il958l ; lAbramowitz fc Stegunl . Il964j ). For m > 0, 
x = —0.1, whereas for m = 0, x = 5.9. The m > model is quite plausible, 
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Fig. 1. Hubble diagram fits for 182 high-confidence SNe la from lRiess et al. I (120071 1. 
The dashed line is a plot of Eq. [3] (flat universe with m = 0) , and the solid line is 
a plot of Eq. [9] (to > 0). The horizontal axis is the redshift z, and the vertical axis 
is the photometrically determined distance modulus. The Hubble constant for each 
curve has been chosen to give the best fit to the observations. 

but the m = model is not; and the estimated H n for the m > model 



is appreciably closer to that of iRiess et al.l (119981 ) than is the estimate for 
m = 0. Indeed, the m > model' s Hubble constant is even clo s er to the 
H Q = 62.3 km s _1 Mpc -1 estimate of iTammann. Sandage fc Reindll ( 20081 ). 



The Einstein-de Sitter (m = 0) and Milne (m > 0) models have power-law 
cosmologies in which the cosmo logical scale factor e volve s as a oc t q . Assuming 
a power law parameterization, Sethi. Dev fc Jair] ( 2005 ) found t he exponent 



that gives the best fit to the first 157 "gold" subset of SNe la data (IRiess et al. 



2004m is q — 1.04^Q Qg. This provides additional support in favor of Milne's 



model (q = 1) over that of Einstein and de Sitter (g = 2/3), but there are 
other cosmologies to consider, notably those involving dark matter and dark 
energy. 



Riess et al.l (120071 ) proffered the ACDM model with cold dark matter and dark 
energy (fi m = 0.29, ^ A = 0.71), fitting their "gold" data set with x 2 = 150, 
which is appreciably less than our Milne fit's x 2 = 179. However, their Milne 
fit's x 2 — 164 is also less than ours; yet they fit the Einstein-de Sitter model, 
with x 2 = 285, which is close to our x 2 = 288. Our calculations differed 
somewhat from theirs, and we could discern no pattern to the offsets, so we 
recalculated their ACDM model fit, and found x 2 — 164, and x = —0.8. Our 
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Table 1 

Qualities of fit between the 182 "gold" supernovae for three cosmologies. All the 
calculations have been independently done for this paper. Where X <C c/H Q = Lh 
(the Hubble length), the m > model is equivalent to a Milne universe. The conven- 
tional cosmological parameters, U m , Ua, and Q#, are associated with the density, 
the cosmological constant, and the curvature; and h is the Hubble parameter. 

Universe $7 m $7 a h x 2 x 

ACDM 0.29 0.71 0.00 0.647 164 -0.8 

Einstein-de Sitter 1.00 0.00 0.00 0.554 288 +5.7 

Milne 0.00 0.00 1.00 0.601 179 -0.1 



fits to the ACDM, Einstein-de Sitter, and Milne models are summarized in 
Table 1. The Hubble parameter, h = if o /(100km s _1 Mpc -1 ), is the only free 
parameter for each fit. The results are independent of the two fitting technique 
we used - marginalization and conventional weighted least-squares. Our h 
estimates are valid only if there are no systematic uncertainties, including 
that of the SN la absolute magnitud e. As an example o f the sensitivity, a 



systematic uncertainty of ±0.03 mag (IHicken et all 120091 ) would cause h to 



be offset by +0.0083. Even so, the systematic uncertainties affect only the h 
estimates, but not the x 2 results. 

If the criterion for the better fit is having the lower x 2 , then the ACDM model 
(x 2 — 164) clearly edges the Milne model (x 2 = 179). However, if we give 
full credence to the observation variances (that is, if we assume that they are 
not exaggerated), then the criterion for the better fit is having the higher 
likelihood. The maximum likelihood is at the peak of the x normal distribu- 
tion, x = 0, so then the Milne model (x = —0.1) edges the ACDM model 
(x = —0.8), but this comparison is for solutions that minimize x 2 rather than 
minimize x. In both cases x < 0, so by varying the solution parameters from 
the least-squares determinations, x will increase from its minimum, and there 
will be non-unique solutions for which x = 0. For example, the only parameter 
for the m > model is h, and the two roots to x(h) = are h = 0.596 and 
h = 0.607. With three independent parameters {h, Q m , fi A ), the ACDM model 
is more complicated because it involves a hypersurface rather than a curve. By 
varying each parameter, one at a time, from the x 2 minimum, six solutions for 
x = can (in principle) be found; but if they are all allowed to vary together, 
then there will be a continuum of solutions. The relative plausibility of these 
two models cannot be definitively settled by statistics. We prefer instead to 
invoke Occam's razor: (1) the ACDM model requires cosmological attractions 
by vast amounts of dark matter that are largely countered by cosmological 
repulsions by vast amounts of dark energy, whereas (2) the m > model 
merely requires that the graviton not be massless. Occam's razor favors the 
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second. Moreover, by not r equiring a A te rm, we avoid conceptual difficulties, 
such as those described by ICarrollI (120011 ). who considers a non-zero A to be 
exceedingly problematic, but who also concedes that denying its existence in 
face of GR and existing data is even more problematic. The massive graviton 
alternative has the additional advantage that it is not sensitive to the gravi- 
ton mass, provided only that m is not identically zero. Finally, if there are no 
massless gravitons, then there is no critical density below which the universe 
is open; the universe is necessarily open. 



4 The graviton mass 



Although the results in Fig. 1 are insensitive to the specific value of m pro- 
vided that m ^ 0, it is in teresting to consider the limits on m that do 
exist. Following IWilll (120011 ) we assume that in the absence of other long 
range forces, the effect of a massive graviton is to replace the Newtonian 
gravitational potential —GM/r of a point mass M by the Yukawa potential 
— (GM /r)e~ r l x . Hence any system of characteristic size R whose behavior is 
correc tly described by Newtonian gravity implies a limi t X > R. Several au- 
thors (IWilj l200ll . [l998l ; iGoldhaber fc Nietol . Il974i : IVisserl Il998h have analyzed 
data for v arious astrophysi c al sys tems to set limits on X. Using solar system 
data from iTalmadge et all (119881 ). Will finds X > 2.8 x 1 12 k m, and from 



galaxy supercluster data one can infer X > 6 x 10 km (IWilll . 120011 . 11998 



Goldhaber fc Nietol . Il974j ; IVisserl . Il998l ). Other lower limits on X are summa- 



rized in IWilll (120011 ). Utilizing the previous arguments we can infer an upper 
limit for X. Suppose that there is an immense volume of uniform density mat- 
ter whose extent is many skin depths (A's). On a point mass that is several 
skin depths inside the outer surface of the volume, the specific force due to the 
Yukawa potential is zero, since the source appears the same in all directions. 
Suppose next that this volume is divided into two parts: one is a sphere of 
radius r < A, and the other part is the remainder of the volume with the 
sphere removed. The sphere is the source of a force directed toward its center 
(which can easily be quantified), so the remainder of the volume is the source 
of a force of an identical magnitude, but directed away from the center of the 
sphere. Hence if there exists a vast uniform density volume, and somewhere 
deep inside it there is a density perturbation that makes the local density 
smaller than the mean density, then the lower density matter will be acceler- 
ated outward. At the outer boundary of the volume matter will, of course, be 
accelerated inwards. The volume will transform toward being a hollow bub- 
ble until the thickness of the higher density bubble matter is just a few skin 
depths. There are some additional considerations: conservation of angular mo- 
mentum places a constraint on how far the inward moving matter can progress 
although it does not, by itself, constrain how far the outward moving matter 
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can go. According to Dolgov. Sazhin fc Zeldovich ( 1990l ). "clusters of galaxies 



form surfaces with a thickness of about 10-20 Mpc surrounding empty regions 
with characteristic size of 100-200 Mpc." With the estimate 

A^5Mpc, (11) 



the gravitational force is significant for galaxy superclusters and also between 
neighboring superclusters that form a "surface," but negligible in the vast 
empty regions between the surfaces. A skin depth A = 5 Mpc = 1.6 x 10 20 km 
corresponds to m = 1.3 x 10~ 30 eV/c 2 ; hence we conjecture that m ~ 10 -30 
eV/c 2 . It is then conceivable that there could be heavier gravitons in addition 
to this lightest graviton, with masses somewhere in the range 

1.3 x 10~ 30 eV/c 2 < m < 7.4 x 10~ 23 eV/c 2 , (12) 



which corresponds to 

1.6 x 10 20 km > X > 2.8 x 10 12 km. (13) 



The Laser Interferometer Space Antenna (LISA) gravitational wave detector 
would not be able to detect the m ~ 10 -30 eV/c 2 mass, but it might be 
sensit ive to the m ore massive gravitons in the range of Eq. [12] if, indeed, any 



exist flWilll . Il998h . 



5 The cosmic microwave background (CMB) 



If m > 0, the dynamics of a baryon fluid at the epoch of recombination 
differ substantially from those of conventional models which tacitly assume 
that the graviton is massless. In Section HI we explained why a massive body 
with dimensions that are large compared with X would not be stable, but 
instead would tend to separate into smaller bodies with dimensions that are 
of the order of a few skin depths or less - that is, to the the scale of galaxy 
superclusters. We hypothesize that this instability is the principal cause of 
CMB power spectrum anisotropics. On recombination, massive bodies with 
breadths larger than, say, D = AX (that is, with radii exceeding two skin 
depths if the bodies are spheres) would divide into lesser sized bodies separated 
from each other by the distance of approximately D. Because the universe has 
expanded, the nominal separation distance of observed radiation peaks from 
the last scattering surface, with redshift z*, is (1 + z±)D; and the nominal 
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angular separation of peaks on the celestial sphere for z+ ^> 1 is (see Eq. [10 
1 + z*)D 2H D 



7 



a>of 



8A 



(14) 



where Lh = c/H a is the Hubble length. Setting h = 0.60 (from Table [3]) and 
X ~ 5 Mpc (Eq. [Till , we calculate 7 ~ 0.008 radians ~ 0.5 ° , whi ch is close to 



what has been observed (ILange et all 1200 ll ; ISpergel et all . 120031 ) . 



6 Summary 



Massive gravitons lead to an alternative relation between the luminosity dis- 
tance [a Q r versus a f from Eq. [10] and the redshift z [Eq. [9] versus Eq. [3]. 
A comparison of Eqs. [9] and [3] in the framework of a cosmology with A = 
favors the massive graviton alternative, as discussed in the text and shown in 
Fig. [TJ Moreover, our graviton mass estimate, m ~ 10~ 30 eV/c 2 , is consistent 
with the structure of galaxy superclusters and with the pattern of the cosmic 
microwave backg round. Although t he redshift data could be attributable to 



20041 ). doing so would conflict with the rule of 



a ACDM model (IRiess et all , 
thumb that "less is better" because of the need for additional suppositions. 



7 Conclusions 



The inverse- square law is the only possible law that leads to a cosmological 
contracting acceleration. Its GR heritage, common acceptance, usage and ap- 
plicability over broad but limited domain lead us to consider it to be the norm. 
It need not be; indeed we might even consider it to be a peculiar exception. 

Any force law that assigns the graviton a mass, m > 0, is a local law because 
the force falls off both geometrically (like the inverse-square law) and exponen- 
tially with distance. In that case, the cosmological principle that "everything 
is the same everywhere" is not even needed. A weaker principle would suffice, 
for example the cosmological legal principle that "every physical law is the 
same everywhere." Gravity pulls the clusters in a supercluster together, but it 
does not pull distinct superclusters together. The gravity-free Milne model is 
underlain by local gravity effects that lead to our interesting universe rather 
than a featureless cloud of expanding matter. Dark matter and dark energy 
are then regarded as gratuitous Procrustean artifacts. 
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A Contrasting the Newtonian and Yukawa Accelerations in a Uni- 
form Density Universe 



We contrast the familiar Newtonian acceleration with the less familiar Yukawa 
acceleration, starting with a model for which there is only a thin spherical shell 
of mass dM at r = a. 



A.l Newtonian Potential Inside the Shell 



The Newtonian potential d^fjsf(r,a) at the origin is d^jv(0,a) = —GdM/a, 
and the Newtonian potential satisfies the homogeneous equation V 2 d^N = 
(in free space), whose solution is a linear combination of d^N = constant and 
d^N oc 1/r. The only such solution that satisfies d^/^(0,a) = —GdM/a is 

dV N (r, a) = -G dM/a, r < a. (A.l) 



The solution interval is closed at the shell because the potential (but not its 
gradient) is continuous there. The potential has no gradient for r < a, so there 
is no Newtonian gravit ational acceler ation inside the shell. This is Newton's 
"iron sphere theorem" (jPeeblesl . Il993l ) . 



A. 2 Yukawa Potential Inside the Shell 



The Yukawa potential at the origin is d^y(0, a) = —GdMexp(—fia)/a. The 
Yukawa potential in free space satisfies the homogeneous equation V 2 <i\I/y — 
fi 2 d^y — 0, whose solutions are linear combinations of exp(±//r)/r. The only 
such combination that equals d^y (0, a) at the origin is d^fy (0, a) sinh(/xr) / (/ir), 
so the Yukawa potential inside and on the shell is 

, T / % ^ ,^ ex P( — I 10 ) sinh(ur) dwAr.a) , „ , . 

d^ Y (r,a) = -GdM— — — ^-2= — il-L- i da, r < a. (A.2) 

a jjr da 



A. 3 Newtonian Potential Outside the Shell 



The solution to V 2 ^^ = that equals d^>N(a,a) at r = a and approaches 
zero as r — > oo is 

cNf N (r, a) = -G dM/r, r > a. (A.3) 
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Outside the shell, the Newtonian gravitational acceleration (the gradient of 
the potential), follows the inverse square law. 



A. 4 Yukawa Potential Outside the Shell 



The solution to V 2 d^/ Y — p 2 d^/ Y = that equals f±(a, a) at r = a and remains 
finite outside the shell is 

, T / \ ^ ,„,sinh(ua) exp(— ur) dw 2 (r,a) , . 

dV Y (r, a) = -G dM ^ PV p ; = — da, r > a. A.4 

pa r aa 

Note that dwi(r,a)/da = dw2(a,r)/da. 

We next integrate over a to derive the Newtonian and Yukawa gravitational 
potentials in a uniform density universe, for which dM = Air pa 2 da. For ^jv, 
this is a trivial exercise. The integral of Eq . IA.3I is —GM/r oc — r 2 , where 
M is the mass of the sphere r < a; this is the familiar Newtonian gravita- 
tional potential. The integral of Eq. IA. 1 1 is —G{Mjj — M)/a, where Mjj is 
the mass of the universe inside the cosmic horizon. The radial gradient of 
the Newtonian potential is GM/r 2 ; at r = a, it proportional to a, so there 
is a gravitational attraction between any two points that is proportional to 
their separation distance. If m = 0, the universe has a uniform contracting 
gravitational acceleration. 



A. 5 Derivation of Eq.{^ 



The integral of Eq. IA.2I is 

= J dw f>«) da = _ 47rGpA2 exp(-H(l + Hsinh(H 
J oa ur 

r 



and the integral of Eq. IA.4I is 

W 2 (r) = ) dW f a) da = —AnGpX 2 -M~^r cosh(pr) - sinh^r)] ^ g) 



Their sum is W x {r) + W 2 {r) = —AnGpX 2 = -3GM A /A, the same as Eq. U If 
m > 0, the universe has no gravitational acceleration. 

If we let p. — > in Eqs. IA.2I and IA.41 then they become Eqs. IA. II and IA.3( 
there is no m = discontinuity for thin shells. However, if we first integrate 
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over all shells, then \I/y is independent of r, whereas oc 1/r; letting /i — > 
after the integration does not change this situation. 



B Derivation of the Metric 



We used Leonard Parker's M athema t ica (IWolframl . ll999l ) notebook, Curvature 
and the Einstein Equation flHartld . |2002|), to determine f(r) in the Eq. Q] 
metric, when the distance formula is Eq. [61 The only change we made to the 
notebook was in the definition of the metric, using these two statements: 

a=aO+adot t 

metric={{-a A 2,0,0,0},{0,-(af[r]) A 2,0,0},{0,0,-(af[r]) A 2Sin^] A 2,0},{0,0,0,c A 2}}. 
Then the G 2 2 {Gee) component of the Einstein tensor turned out to be 



G 



22 



f(r) f"{r) - (a/cf 



(B.l) 



Eq. [7| is an evident solution to G22 = 0, so we substituted it back into the 
Mathematica notebook by defining 

f[r] = (c/adot) Sinh[adot r/c], 

and all components of the Einstein tensor turned out to be zero. QED. 
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